Optical conductivity of Mn doped GaAs. 
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The optical conductivity in the III-V diluted magnetic semiconductor compound (Ga,Mn)As is 
investigated theoretically. A direct comparison between calculations and available experimental 
data is provided. We demonstrate that our model study is able to reproduce both qualitatively and 
quantitatively the measurements performed by Burch et al. and Singley et al.. It is found that the 
increase of the carrier density up to one hole per Mn leads to a red shift of the broad conductivity 
peak located at approximately 200 meV in optimally annealed samples. Our study demonstrates 
that the non perturbative treatment (beyond the valence band picture) is crucial to capture these 
features. Otherwise a blueshift and an incorrect amplitude would result. We have calculated the 
Drude weight (order parameter) and have established the metal-insulator phase diagram. It is shown 
that (i) Mn doped GaAs is indeed close to the metal-insulator transition and (ii) in both 5% and 7% 
doped samples, 20% of the carriers only are delocalized. Beyond the mobility edge, we have found 
that the optical mass is m op t ~ 2 m e . Interesting new features for overdoped samples have been 
found. The overdoped regime could be experimentally realized by Zn codoping. Detailed discussions 
and carefull analysis are provided. 



The possibility to use for spintronics purpose the spin 
degree of freedom in magnetic semiconductors has lead 
to a tremendous experimental and theoretical work. One 
of the main goals is to search for optimal candidates that 
exhibit high Curie temperatures (beyond room tempera- 
ture). Although it is well accepted that the presence of 
holes in III-V diluted magnetic semiconductors is respon- 
sible for the induced ferromagnetism, the particular case 
of (Ga,Mn)As remains still controversial concerning the 
nature of the states at the Fermi level. There are in fact 
several important points that arc still controversial in the 
literature. Is the Fermi level inside the weakly perturbed 
valence band (VB scenario) or does (Ga,Mn)As exhibit 
a well defined impurity band ? Are the holes extended 
or localized ? Or is the situation similar to that of Mn 
doped II- VI as (Zn,Mn)Te ? For those latter compounds 
the Valence Band picture is appropriate. It seems that 
this could be at the origin of some of the controversies. 
In our view the main question should be: In the case of 
(Ga,Mn)As, is a perturbative treatment (Valence Band 
picture) appropriate to explain the physics ? Or equiv- 
alently, is it necessary to treat disorder effects properly? 
From ab-initio based studies, one finds that (Ga,Mn)As 
clearly exhibits a preformed impurity bandi~— that is re- 
sponsible for (i) the non RKKY nature of the magnetic 
couplings and (ii) explains the measured high Curie tem- 
peratures (with respect to II- VI compounds). Addition- 
ally, it has been clearly demonstrated that RKKY cou- 
plings obtained within the VB picture can not explain 
the magnetism in III-V compounds^. Indeed, the frus- 
tration effects would lead to either a spin glass phase or 
to extremely small Curie temperatures. On the other 
hand, the non perturbative treatment of the appropri- 



ate model Hamiltonien captures the essential physics and 
lead to non RKKY couplings^, in agreement with ab ini- 
tio studies. Note also that, below 1% the impurity band 
is completely splittcd from the valence band. Early re- 
sistivity measurements have clearly demonstrated that 
(Ga,Mn)As is close to the metal-insulator phase transi- 
tion. Indeed, measurements in as grown samples often 
show an insulator behaviour at low temperature and a 
metallic one after annealing^—. Of course, the VB sce- 
nario is unable to capture this physics, namely the metal- 
insulator transition. It has been shown, by Singley et al. 
and Burch et al. that the measured optical conductivity 
is also inconsistent with the perturbative Valence Band 
scenario^—. The aim of this paper is to provide a non 
perturbative theoretical study able to reproduce and ex- 
plain the experimental observations. 

Ab-initio based studies have provided the most reli- 
able tool, to allow for quantitative studies of the mag- 
netic properties without any adjustable parameters. As 
an example, it was possible to study in great details and 
quantitatively the magnetic properties of Gai_a;Mn x As, 
both in the presence and absence of native (compensat- 
ing) defects^ - — . However, dynamical transport study 
is not a simple task within the ab-initio approach. The 
study of relevant minimal model is thus necessary. The 
model approach is especially suitable to understand the 
influence of a particular physical parameter and may pro- 
vide a support for finding new potential spintronic can- 
didates. The V-J model treated non perturbatively was 
shown to capture qualitatively the magnetic properties 
of a whole family of III-V materials^. The present the- 
oretical study is based on such a model. The one band 



2 



V-J Hamiltonian reads, 

H = - ^ tijcl a Cj a + };PiJjSj ■ S 4 + 2J ^PiC^Cicr (1) 

In the first term, ty=t if i and j are nearest neighbors, 
otherwise ty=0. The random variable ^ = 1 if the site is 
occupied by an impurity, otherwise it is 0. denotes the 
localized impurity spin at site i (\Si\ = 5/2) and the 
spin 1/2 operator of the carrier. is the p-d coupling 
(J) between itinerant carrier (p-states) and localized Mn 
spin . The last term results from the substitution of Ga 3+ 
by Mn 2+ : ei =V if the site is occupied by Mn, otherwise 
it is 0. This additional crucial term allows to adjust the 
position of the hybridized p-d states with respect to the 
top of the VB^^. In the following, x denotes the Mn 
concentration and p the hole density. In the absence of 
compensating defects (As anti-sites or Mn interstitials) 
p = x. In order to perform reliable calculations, the local 
coupling between impurity and carrier is treated in a non 
perturbative way and the disorder and dilution effects ex- 
actly (no effective medium approach) . This treatment al- 
lows for the localization of the itinerant carriers which is 
expected to strongly affect both magnetic and transport 
properties. The procedure is as follow. For a given impu- 
rity concentration x and a given disorder configuration, 
the calculations are done by exact diagonalization of the 
Hamiltonian (TTJ) in both spin sectors. This provides the 
set of eigenvalues and eigenstates {E°' c , \^' c )} where 
<7 =t,4- an d r = 1,2, N , that are used to evaluate 
the optical conductivity <t{uj) and the Mn-Mn magnetic 
couplings. The superscript c denotes the configuration of 
disorder, it will be omitted in the following. Note that, 
the calculations are performed at T= K assuming that 
the impurity spins are aligned. 




FIG. 1: (Color online.) Illustration for uncompensated sam- 
ple. In our model (left), one impurity brings one state above 
the top of the valence band. In reality (right), each Mn pro- 
vides rii hybridized p-d states (nj = 3) and 1 hole/Mn. To be 
consistent, in the model the hole concentration is thus set to 
p = x/m. 

For simplicity, our calculations are performed on a sim- 
ple cubic lattice instead of the fee lattice of Ga in the real 
GaAs compound. Thus, instead of 4 atoms/unit cell we 
have only one. Hence, the value of the lattice spacing 
that will be used is a — -jf^ = 3.55 10~ 10 m where 
ao = 5.65 10~ 10 m is the lattice parameter of the Zinc- 



Blende GaAs material. Let us, shortly justify this choice. 
In these diluted compounds the carrier densities are very 
low, thus it is natural to expect the structure of the lattice 
to play no relevant role. As will be seen, this simplifica- 
tion will not have a significant effect on the results. This 
is also valid when the carrier are localized, as far as the 
localization length remains much larger than the lattice 
spacing. 

Let us now discuss how the model parameters are fixed. 
First, the hopping integral t is chosen according to the 
bandwith W of GaAs host. Thus, we use t « 0.7 eV 
(W=12 t). JS is set to 3 eV since it is now well ac- 
cepted that in (Ga,Mn)As J«1.2 eViZ»!£. Note that J 
is also of the order 1 eV for other materials including 
II- VI compounds ((Zn,Mn)Te, (Cd,Mn)Te,..)^^. The 
last parameter V (on-site energy) is as said above the 
crucial one. It actually explains the different nature of 
the Mn-Mn couplings and Curie temperatures values in 
Mn doped compounds such as (Zn,Mn)Te, (Ga,Mn)As 
and (Ga,Mn)N. The on-site energy V contains two con- 
tributions: (i) an electrostatic one due to the difference 
of charge between the host and substituted cation and 
(ii) another resulting from the Schricffcr- Wolff transfor- 
mation of the Anderson Hamiltonian 2 ^. Indeed, in ad- 
dition to ^iPiJiSi ■ Si term, this transformation leads 
to a spin independent scattering contribution which de- 
pends on the position of the d-levels with respect to 
the top of the valence band in the host material. V 
is set in order to reproduce the energy of the accep- 
tor level E 6 =110 meV in GaAs host^-22. This leads 
to V=1.25 cV=1.8 t. Note that with this set of pa- 
rameters, we were able to reproduce the whole varia- 
tion of the Zeeman splitting as a function of x obtained 
within first principle studie d 29 ' 30 . This splitting is defined 
by A(x,V) = i4 ax - wh ere ££ ax is the largest 

eigenvalue in the cr-sector—. For example, from ab-initio 
calculations^ -, A ab " initio = 0.65 eV for x = 0.05, we 
obtain a value of 0.66 eV. Thus, from now on, the set 
of three parameters (t, J, V) is fixed and the calcula- 
tions that will follow are done without any adjustable 
additional parameter. The calculations are performed for 
both 5% and 7% doped compounds in the ferromagnetic 
phase. Experimental and theoretical studies have shown 
that below 1% no long range ferromagnetic order is possi- 
ble (below the percolation threshold). In the very dilute 
regime (below 1%), the optical conductivity was recently 
studied in refi^i within a very different theoretical ap- 
proach. They perform some variational calculations on 
Mn-Mn molecule and extract the model parameters that 
are used to perform their calculations. 

Let us now explain in more details our theoret- 
ical method. The total optical conductivity reads 
c(<j) = X)<T Cr o-(^) where <J a {u)) is given by the follow- 
ing expression 3 ^—, 

a tr (u>) = D (r S(uj) + a r ^(u ) ) (2) 

D a is the Drudc weight (w = contribution) in the 
cr-sector and a r a eg (u>) is the regular part of the optical 
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conductivity. Note that the Drude weight is the order 
parameter for the metal-insulator phase transition: it is 
zero in the insulating phase and finite in the metallic 
one. In our calculations, we have used periodic bound- 
ary conditions in order to separate the finite frequency 
contribution from the dc part. In that case, the Kubo 
expression provides <7 r a eg (uj) only. The Drude weight D is 
then obtained from the following sum-rule, 

(K® ) is the average hole Kinetic energy in the x direction 
and CT-sector, N=L 3 is the total number of sites, and 
(7 = t- = 21540 fi-^cm- 1 . 




FIG. 2: (Color online.) Illustration for the procedure of cal- 
culation of the magnetic couplings between two Mn impuri- 
ties located at site i and j. The left part corresponds to the 
real compound in which each Mn 2+ provides 3 p-d hybridized 
states and one hole (this picture describes 3 Mn impurities 
in GaAs). The right side corresponds to the one band V-J 
model. The c^J are the matrix elements that enter the ex- 
pression of the exchange given in eq. ©. 

To allow for a direct comparison between available ex- 
perimental data and theory, we have to take into account 
the following point: the substitution of Ga 3+ by Mn 2+ 
provides n/ = 3 p-d states near the top of the valence 
band (see FigfT] for an illustration). These states are 
known as the Dangling Bond Hybrid states 3 -^. As men- 
tioned before, in the case of a single Mn impurity, the 
acceptor level energy is Ef, rj 0.11 eV. At finite con- 
centration, an impurity band is formed and merges into 
the VB when x becomes large enough x > 0.010. Be- 
cause Mn 2+ provides 1 hole per Mn, the resulting p-d 
hybridized states impurity band is thus 1/3 filled. On 
the other hand, in the one band V-J model, each impu- 
rity provides a single state. Thus, to be consistent and 
coherent, if the (Ga,Mn)As compound contains a concen- 
tration x of Mn and p of holes the theoretical calculations 
should be done as follows. (1) The one band calculations 
are performed for the hole density p = — and (2) the 
theoretical optical conductivity that is compared to the 
experimental data is, 

<r(u,p) = nia(u,p/ni) (4) 

where d(w,p/n{) is that calculated within our one band 



model. Similarly, the Curie temperature will be per- 
formed using the Mn-Mn magnetic couplings, 

Jij{x,p) =niJi,j(x,p = p/ni) (5) 

where, Ji^{x,p/n{) are those calculated within the one 
band V-J model. To be more specific these couplings are 
given by the following generalized susceptibility, 

(6) 

where the one particle retarded Green's functions 
Gf„-(w) = (i(j\ — i — h'cr). In the present case, the lo- 
cal exchange splitting energy is a constant S, =J p( jS. An 
illustration of the procedure that provides the Mn-Mn 
exchanges is shown in Fig|2] Note that, this picture ap- 
plies also for the optical conductivity calculation. In this 
case, the matrix elements c|J are replaced by those in- 
volving the current operator (see below). The previous 
equation can be rewritten as follows, 

r,r' r T 1 

where n^ is the occupation number associated to the 
state \%Z). The coefficients df = A], r .A^ r ,, where 

Afj, r = The matrix elements c[J con- 

tains the relevant informations concerning the nature of 
the hole states, thus the nature of the magnetic cou- 
plings. In the VB picture r (resp. r') are replaced by the 

momentum k (resp. k') and c-J = le'f*:-' )'( r .-rj) 
leading to the standard RKKY couplings. The diag- 
onalization of the resulting effective dilute Heisenberg 
Hamiltonian Hneis = ^ijPiPjJi,j{ x iP)Si'Sj will provide 
the Curie temperature Tc{x,p). This procedure will be 
done within the self-consistent local RPA (SC-LRPA)^. 
The self-consistent local RPA is a semi-analytical Green's 
function based method which is in particular able to cap- 
ture quantitatively both the physics of localization and 
percolation. The disorder is treated in the real space 
(no effective medium) and the thermal/transverse fluc- 
tuations beyond standard Mean Field within the RPA. 
Comparisons with Monte Carlo treatment have demon- 
strated several times that the SC-LRPA is a reliable and 
very accurate tool (see for exampl o 14 ' 36 ). 

The regular part of the optical conductivity 
a^ 9 (w , p j ni) , is given by the Kubo expression which 
reads, 

K e9 (",P/ni) = f E« - <') E fl' E J ^ E r + *v)(8) 

The matrix element AZ r , = |(#£|j£|*p)| 2 , where j£ = 

—itJ2ij( c lja c i+x,jcr — he) is the x component of the cur- 
rent operator in the cr-sector. 
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FIG. 3: (Color online.) Regular part of the optical conduc- 
tivity cr(w) as a function of uj/t for x — 0.05: (left) for 'under- 
doped' samples (electron doped) , e.g. from low carrier con- 
centration to p/x = 1 (well annealed sample) , (right) 'over- 
doped' samples (hole doped). The parameters are JS=4.3 t 
and V=1.8 t ((Ga,Mn)As set). The optical conductivity is 
given in units of ctq (see text). 




p/x < 1 1 < p/x < 2 2 < p/x < 3 



FIG. 4: (Color online.) The top of the figure is a sketch of 
the density of states for Mn doped GaAs: a compound with 
contains typically 5% of Mn. The zoom of the preformed 
impurity band is also shown. There are three relevant regions 
which provide the dominant transition matrix elements in the 
Kubo formula of the optical conductivity (see eq.((8|). There 
are 3 relevant regions. Region (1) corresponds to the localized 
states (the mobility edge Ec is indicated) and both Regions 
(2) and (3) to disordered extended states, (a) For p/x < 1 
the dominant transition matrix element is A 12 between region 
(2) and region (1); (b) For 1 < p/x < 2, the matrix elements 
A12 reduce and A23 (between region 2 and 3) are now allowed 
and start to dominate; (c) for 2 < p/x < 3, A12 contributions 
are not possible anymore and the dominating process is A23. 
Note that the width of the vertical arrows indicates the weight 
of the processes. 



In Figj3] we have plotted, for x = 0.05 and various 
carrier densities the regular part of the optical conduc- 
tivity as a function of frequency. The calculations are 
realized on a 20 3 sites system. A comparison with those 
done on both 16 3 and 24 3 systems has shown that the fi- 
nite size effects are negligible. A systematic average over 
at least 200 configurations of disorder has been done. 
As observed experimentally, a (to) exhibits a clear broad 
peak for all carrier concentrations. As the hole density 
increases from the lowest concentration up to optimal 
doping p/x = 1 (left panel of FigJ3]), we observe a clear 
redshift of the peak position. This feature is in agreement 
with the detailed experimental study performed by Burch 
et al. (see left panel of Fig. 3 in refJ^). Moreover, this is 
in contrast to the valence band picture which predicts a 
blueshift 3 ^ (see also right panel of Fig. 3 in ref>^). Note 
that, a blue-shift was reported in recent measurements 3 ^ 
(see Fig.l (c) in the preprint), but no clear explanation 
was given to explain the origin of the conflict between 
earlier studies and their recent results. The data shown 
in Fig.l (c) in ref<2& result from a fit procedure which 
depends on many parameters. Thus, it would be of great 
interest and very useful if the authors could show the 
measured data, this would facilitate the comparison with 
earlier studies. On the other hand, we have noticed that 
the annealed samples appear to be well annealed. In- 
deed, the measured Curie temperatures Tc (shown in 
Table I of refi^£) agree very well with (i) the first princi- 
ple based calculations^ and with those obtained in the 
present study (see below). 

From the left panel of Fig(31 we observe that the peak 
is located at 0.2 t for the uncompensated sample (p = x) 
and approximately 0.35 t for the lowest concentration 
which corresponds to the presence of about 25% of com- 
pensating native defects as As anti-sites (double donors 
of electrons) for example. A direct quantitative com- 
parison with the experimental data will be given in the 
next section. To our knowledge, no experimental data 
are available for larger hole density (beyond 1 hole/Mn), 
nonetheless we have also analyzed the effect of co-doping 
(Ga,Mn)As with a hole donor. This can be achieved ex- 
perimentally by the substitution of Ga 3+ by Zn 2+ for ex- 
ample (Zn 2+ introduces one hole/Zn). We observe in the 
right panel of Fig J3] that adding holes leads to an increase 
of the conductivity and now to a blueshift of the peak to 
much higher energy. For example for p/x > 2.1, the peak 
is now located at approximately l.Oi. Such a hole con- 
centration could be obtained in Gao.9Mno.05Zno.05As for 
example. It would be interesting to analyze experimen- 
tally the variation of <j(uj) as a function of y in the se- 
ries of compounds Gao.gs-yMno.osZn^As after annealing. 
Let us now discuss in more details the origin of both the 
redshift and blueshift obtained. For different hole density 
concentrations, we have been able to isolate the processes 
which dominate in the optical conductivity response. In 
other words, we have identified which matrix elements 
A a rr , are responsible for the peak in <r re9 (uj,p). We could 
identify 3 relevant regions denoted (1),(2) and (3) de- 
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picted in Fig|4] The origin of the red-shift is explained 
in the cartoon (a): below p/x =1, the only possible tran- 
sitions are between region (1) and (2) (transitions are 
only allowed between occupied and empty states). The 
transitions between (1) and (3) are much smaller. As the 
hole density increases up to p/x = 1, the Fermi energy 
moves down and the energy W12 associated to A12 pro- 
cesses reduces from 0.4 t to 0.2 t. This is the explanation 
of the red-shift. In the very dilute limit (typically below 
1%), this corresponds to the process between the well 
defined (splitted) impurity band and the top of the va- 
lence band as discussed in refJ^. For 1 < p/x < 2 (case 
(b)) the Fermi level lies now in the region (2). The two 
relevant matrix elements are of type A12 and A23. As 
the hole density increases, the relative weight of A23 in- 
creases whilst A12 reduces. This can be seen in the right 
panel of Figj3]for p/x = 1.35 and 1.95, where a shoulder 
coming from the A23 processes is clearly visible and gets 
larger as the hole density increases. The typical energy 
associated to A23 processes is of the order W23 ~ t. As 
we further increase the hole density (case(c)), A12 pro- 
cesses are now forbidden (regions (1) and (2) are empty 
or filled with holes) and the only possible transitions are 
between the regions (2) and (3). Thus, this explains the 
peak located at j w 1 for p/x > 2 and the origin of the 
blue-shift. 

We now proceed further and discuss the variation of 
the metal-insulator order parameter, the Drude weight, 
as a function of the carrier concentration. Note that, the 
dc conductivity can be obtained from <jdc = Dt, where 
t = l e /vp, l e is the mean free path. In the top panel 
of Fig 13 the averaged Drude weight Dl is plotted as a 
function of the carrier density per Mn for various system 
sizes. The Mn concentration is set to x = 0.05 and L 
varies from L—12 to 20. For the smallest system an av- 
erage over 1000 configurations of disorder has been done 
and for the largest 200 configurations were used. Af- 
ter performing an extrapolation in the thermodynamic 
limit (see the bottom panel of figE]), we have obtained 
the metal-insulator phase diagram plotted in the same 
figure. As seen, below a critical concentration (mobility 
edge) the Drude weight D is zero: the system is an in- 
sulator. Above this threshold it scales linearly with the 
delocalized (extended) hole concentration, p cx — (p—p c ). 
The transition point is located at (£) c ps 0.80 for both 
Mn concentrations x = 5% and x = 7.3%. Thus our 
present study predicts that the critical carrier density is 
p c ps 0.0375 for the 5% doped sample and p c ps 0.055 
for x = 0.073. This result is important since it indi- 
cates that approximately 20% of the carriers only con- 
tribute to the Drude weight, most of them (80%) are 
localized. Beyond the critical concentration, the slope 
of the Dude weight D provides the value of the opti- 
cal masses (for localized carriers m opt = 00). Indeed, 

in this regime, one can write D ^ p ' x ^ = Pcx6 , We find 

m op t(x = 0.073) ps m opt (x = 0.05) ~ 4 mo = f , where 
mo ~ 0.5 m e is the hole effective mass in GaAs host. All 
these results are consistent with the experimental obscr- 
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FIG. 5: (Color online.) (top panel) Drude weight Dl (in 
units of t) as a function of the carrier density per Mn, p/x, 
calculated for various systems sizes (L=12,14,....20), with x = 
0.05. (bottom panel) Extrapolation in the thermodynamic 
limit: Doo, the shaded region corresponds to the insulating 
phase, two different Mn concentrations are shown x = 0.05 
and 0.073. The material parameters are respectively JS = 
4.3t and V = 1M 



vation that (Ga,Mn)As is close to the metal insulator 
phase transition and why as grown insulating samples 
become metallic after annealing. Note that from their 
experimental data Singley et al. have estimated m opt to 
be within the range 0.7m e < m opt < 15m £ r^. 

The natural question which arises now is the following: 
Is our theoretical model also able to explain the experi- 
mental data quantitatively? Let us then compare directly 
our calculations to the experimental data of the optical 
conductivity measured in annealed samples^. We will 
see that our model theory is able to reproduce the optical 
conductivity data and the Curie temperatures simultane- 
ously. In Fig|Bl we have plotted the optical conductivity 
in absolute units (f2-cm -1 ) as a function of the energy 
(in eV). The experimental data for the annealed samples 
are taken from refJ^. The calculations were done assum- 
ing one hole per Mn (perfectly annealed samples) , taking 
into account, as discussed before, the fact that each Mn 
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FIG. 6: (Color online.) Calculated and measured optical con- 
ductivity a{u) in Gai-zMn^As (in units S7 _1 -cm _1 ) as a func- 
tion of oj/t for the optimal doping (lhole/Mn) and for both 
Mn concentration x = 0.05 and x — 0.073. 

impurity brings three hybridized p-d states. We observe 
a very good and surprising quantitative agreement for 
both concentrations of Mn. It is even excellent in the 
case of the 5% doped sample. The peak position is lo- 
cated at 200 meV for x = 0.05 in perfect accordance with 
the measurements and the overall shape (peak structure, 
width) is also very well reproduced. Furthermore, using 
the calculated couplings Ji,j(x,p) in the effective diluted 
Heisenbcrg Hamiltonian (see above), we have calculated 
within SC-LRPA the Curie temperature for both cases. 
We have found that Tc is respectively 120 K and 150 K 
for x = 0.05 and 0.07. These values agree very well with 
those measured in these annealed compounds, namely 
120 and 140 K respectively-^. Note also that our cal- 
culated values agree very well with those obtained us- 
ing the ab-initio couplings^. It is interesting to remark 
that if instead of a = -^jj we would have used do, then 
the calculated optical conductivity would have been 60% 



smaller. The small deviation between theory and experi- 
ment observed in Figj6]for x = 0.073 may be explained as 
follows. We have used for the comparison in both cases, 
the total measured optical conductivity (Drude part in- 
cluded). As seen in Fig.l of refJ^ a small shoulder is 
seen at 30 meV. The Drude contribution is larger in the 
x = 0.073 compound than in the 5% doped one and thus 
this may explain why below the peak, the measured op- 
tical response is a bit larger than the calculated one. 

To conclude, we have calculated in a non perturbative 
way the optical conductivity, the phase diagram and 
Drude weight as a function of the hole density in GaAs 
samples doped with 5 % and 7.3% of Mn. In agreement 
with experimental measurements^, we have obtained a 
redshift of the optical conductivity peak. We were also 
able to reproduce quantitatively the measured optical 
conductivity in well annealed samples. Additionally the 
calculated Curie temperatures are in excellent agreement 
with those obtained using ab initio couplings and with 
the experimental values. This success shows that our 
minimal model is able to capture not only qualitatively 
but also quantitatively both transport and magnetic 
properties in diluted magnetic semiconductors. 
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